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The thermodynamics of layered antiferro- and ferromagnets with a weak inter- 
layer coupling and/or easy-axis anisotropy is investigated. A crossover from an 
isotropic 2D-like to 3D Heisenberg (or 2D Ising) regime is discussed within the 
renormalization group (RG) analysis. Analytical results for the the (sublattice) 
magnetization and the ordering temperature are derived are obtained in different 
regimes. Numerical calculations on the base of the equations obtained yield a good 
agreement with experimental data on La2Cu04 and layered perovskites. Corre- 
sponding results for the Kosterlitz-Thouless and Curie (Neel) temperatures in the 
case of the easy-plane anisotropy are derived. 
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The problem of layered magnetic systems is of interest both from theoretical and 
practical point of view Here belong, e.g., quasi-two-dimensional (quasi-2D) per- 
ovskites, ferromagnetic monolayers and ultrathin films. We consider the Heisenberg model 
with small parameters of the interlayer coupling a = J'/ J and easy-axis anisotropy 
T] = 1 — / J^^y (J is the in-plane exchange parameter). This case permits a regular 
consideration since the magnetic transition temperature is small, Tm <^ |<^|'S'^, which is 
reminiscent of weak itinerant magnets. 

At not too low temperatures T the standard spin-wave theory (SWT) is insufficient 
to describe correctly thermodynamics of such systems. Somewhat better results can be 
obtained within the self-consistent theory (SSWT) p, 0, which takes into account the 
temperature renormalization of a and r]. However, the values of the ordering temperature 
in SSWT are still too high in comparison with experimental ones, and the critical behavior 
is quite incorrect. 

To improve radically SSWT, the summation of leading contributions in all orders of 
perturbation theory should be performed. To this end we use the renormalization group 
(RG) approach similar to that of Ref. @|. This approach is valid outside the critical 
region which is very narrow for layered systems. The same results can be obtained by 
direct summation of RPA-type corrections to spin-wave interaction vertex. 

The result for the relative (sublattice) magnetization ct,. = S/S(T = 0) reads [0, || 
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Af], other quantities are given in the Table 1: 
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where c = y/SlJljS is the spin-wave velocity, 7 ~ 1 + 0.078/5" is the renormahzation 
parameter for intralayer couphng, Zli = Zl2 = Zl^ = 1 — T/87r| Jj^^. The inequahty 
r A should to be satisfied for the vahdity of ([1|). Note that the classical regime 
is realized only for very large S. The quantities Jt and in (|l]) are the temperature- 
renormalized values of interlayer coupling and anisotropy parameters, and for T <^ AnpsWr 
(i.e. beyond the critical region, see below) 

fr/fr = {ar/arf = wl (2) 

Thus the anisotropy and interlayer coupling are strongly renormalized with temperature 
which should be taken into account when treating the experimental data. In the quantum 
regime, the parameters and fr are ground-state (quantum-renormalized) anisotropy 
and interlayer coupling (for ferromagnetic case the ground-state renormalizations are ab- 
sent). In the classical regime, fr and ar are 'lattice-renormalized' parameters of anisotropy 
and interlayer coupling, i.e. the corresponding parameters of the continuum model with 
the same thermodynamic properties as the original lattice model. 
Three temperature regimes can be distinguished. 

(i) low temperatures, T ^ Tm ~ 2'k\J\S^/ ln(r/A). The analysis of non-uniform suscepti- 
bility shows that the excitations in the whole Brillouin zone have spin-wave nature. Only 
first term in the square brackets in (|l|) is to be taken into account and the magnetization 
also demonstrates the spin- wave behavior. 

(ii) intermediate temperatures, (S'/S')/ln(r/A) ^ T/27i\J\S'^ ^ S/S (T is of the same 
order as Tm)- Close to the center of the Brillouin zone the excitations still have the 
spin-wave nature, while for large enough momenta they have non-spin-wave character. 
Only in-plane (two-dimensional) fiuctuations are important in this regime. All the terms 
in (|I]), except for the last, are important, which leads to significant modification of the 
dependence ^^(T) in comparison with SWT. 

(iii) critical region, T/ {27t\J\S'^) ^ S/S (1 — T/Tm <^ 1). In this regime the spin-wave 
excitations are present only for momenta ^ A (hydrodynamic region) whereas at all 
other q excitations have non-spin-wave character. The thermodynamics in this case is 
determined by 3D (or Ising-like) fiuctuations. The contribution of $a is of the same order 
as other terms in the square brackets of (^ and the RG approach is unable to describe 
the thermodynamics in this regime. The values of critical exponents can be corrected in 
comparison with SSWT with the use of the expansion. For the isotropic quasi-2D 
antiferromagnet we obtain H] 
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with Aq ~ 0.9635 and /^a = (1 - %/'k'^N)/2 ~ 0.36. 

Up to some (unknown) constant C{f /a) we have for the transition temperature the 
result 

Tm 
Airp 

Note that all the logarithmic terms are included in (^) and C gives only a small contribu- 
tion to this result. For the quantum antiferromagnets the value of C{oo) can be calculated 
by the 1/N expansion 0, C{oo) ~ —0.0660. The value of C(0) for the same case can be 
deduced from experimental data on layered magnetic compounds 0, C(0) ~ —0.7. 
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For practical purposes, simple interpolation expressions for the functions which 
permit to describe the crossover temperature region, are useful. We obtain at a; < 1 
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Numerical calculations with the use of the equations obtained yield a good agreement 
with experimental data on the layered perovskites (see Figs. 1-2 and the Table 2), and 
the Monte-Carlo results for the anisotropic classical systems. 
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In the easy-plane case [rj < 0) finite-temperature magnetic transition is absent at 
a = 0. At the same time, the Kosterlitz-Thouless transition where unbinding of topological 
excitations (vortex pairs) occurs. For small anisotropy value the temperature of this 
transition is small, and to leading logarithmic accuracy M 
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which is similar to the result for Tm in the easy-axis case. As well as for magnets with 
small easy-axis anisotropy, topological excitations are important only close to Tkt- Using 
the renormalization group approach, the result similar to (HI) can be obtained H 
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For the magnetic ordering temperature in the presence of interlayer coupling we obtain 
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The comparison of our results with the experimental data is presented in the Table 3 
(the parameters C ~ —1.0 and A = 3.5 are fitted for the first compound). 



Compound 
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11.4 
5.5; 6.25 
5.5; 6.25 


35.3 
17.4; 18.7 
18^20 


45.0 
23.0; 24.3 
23.0; 24.5 



One can see that the two- loop corrections improve radically the result (|^). 
Figure captions 

Fig.l. The theoretical temperature dependences of the relative sublattice magnetiza- 
tion ar from the spin- wave theory (SWT), SSWT, Tyablikov theory (TT), RG approach 
(0), and the experimental points for La2Cu04 0. The RG' curve corresponds to inclusion 
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of the function ^^^{t/a) given by (^. The curve is the critical behavior predicted 
by 1/iV-expansion @. The resuh of 1/A^-expansion at intermediate temperatures 
practically coincides with RG'. 

Fig. 2. Temperature dependence of the relative sublattice magnetization a(T) of 
K2NiF4 in the SWT, SSWT, RG approaches and 1/A^-expansion for 0{N) model as 
compared with the experimental data [|1^ (circles). The RG' curve corresponds to inclu- 
sion of the function ^^^(t/a) given by (^. Short-dashed line is the extrapolation of the 
1/A^-expansion result to the critical region. 
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